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O " ABSTRACT 

<N 

Context. The nearby young stellar association e Cha has an estimated age of 3-5 Myr, making it an ideal laboratory to study the disk 
O 1 dissipation process and provide empirical constraints on the timescale of planet formation. 

Aims. We wish to complement existing optical and near-infrared data of the e Cha association, which provide the stellar properties of 
its members, with mid-infrared data that probe the presence, geometry, and mineralogical composition of protoplanetary disks around 
individual stars. 

(N- Methods. We combine the available literature data with our Spitzer IRS spectroscopy and VLT/VISIR imaging data. We use proper 

motions to refine the membership of e Cha. Masses and ages of individual stars are estimated by fitting model atmospheres to the 
optical and near-infrared photometry, followed by placement in the HR-diagram. The Spitzer IRS spectra are analyzed using the 
two-layer temperature distribution spectral decomposition method. 

Results. Two stars previously identified as members, CXOU J120152.8 and 2MASS J12074597, have proper motions that are very 
different from those of the other stars. But other observations suggest that the two stars are still young and thus might still be 
related to e Cha. HD 104237C is the lowest mass member of e Cha with an estimated mass of ~ 13— 15 Jupiter masses. The very low 
mass stars USNO-B 120144.7 and 2MASS J12005517 show globally depleted spectral energy distributions, pointing at strong dust 
settling. 2MASS J12014343 may have a disk with a very specific inclination, where the central star is effectively screened by the 
cold outer parts of a flared disk, but the lOyum radiation of the warm inner disk can still reach us. We find that the disks in sparse 
stellar associations are dissipated more slowly than those in denser (cluster) environments. We detect C9H2 rovibrational band around 
13.7 /jm on the IRS spectrum of USNO-B 120 144.7. We find strong signatures of grain growth and crystallization in all eCha members 
with 10 fim features detected in their IRS spectra. We combine the dust properties derived in the e Cha sample with those found using 
identical or similar methods in the MBM 12, Coronet, rj Cha associations, and in the cores-to-disks legacy program. We find that disks 
around low-mass young stars show a negative radial gradient in the mass-averaged grain size and mass fraction of crystalline silicates. 
A positive correlation exists between the mass-averaged grain sizes of amorphous silicates and the accretion rates if the latter is above 
~10~ 9 M yr _1 , possibly indicating that those disks are sufficiently turbulent to prevent grains of several microns in size to sink into 
the disk interior. 

OO 

l/") 1 Key words, open clusters and associations: eCha - stars: pre-main sequence - planetary systems: protoplanetary disks 

Ov 
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CN| 1. Introduction cies in the MBM 12 and ^Cha associations suggest that disk 

■ dissipation proceeds more slowly in these sparse P MS associ- 

•' ■ The observational characterization of the structure and evolution ations compared to denser environments (see, e.g., iFang et al.1 

• H.° f circumstellar disks is key to our understanding of the disk [20I2I) . MBM 12 has a di sk frequency of ^67% at an age of 

K> . dissipation and planet formation processes. Nearby pre-main se- 2 Myr (iMeeus et al.ll2009l while in the 5- 10 M yr old 77 Cha as- 

5^ ; quence (PMS) associations at ~50-300pc are well suited for de- S o C iat ion, 40-50% of the sta rs still retain a disk ( Mamaie ketall 

Cd . tailed investigations of young stars and their disks, since their |l999l iMegealhiTall 120051 Isicilia-Aeuilar et al.1 120091) . In the 

members can be observed with high signal-to-noise ratio (S/N) latter clusterj binarity was sh own to play an impor tant role in 

throughout the electromagnetic spectrum. Several PMS associa- prot oplanetary disk evolution (iBouwman et al.ll2006h . 

tions have been studied with the Spitzer Space Telescope, e.g. Silicates are an important dust component in protoplan- 

MBM12, eCha, and nCha (Meeus et al. 2009; MegeathetalJ e tary disks. The most convincing evidence that silicates are 

|2005t I Bouwman etal. | | 2006t [Sicilia-Aguilar et al. 2009) using present in me p rot0 pi ane tary disks is the strong "10/an fea- 

the infrared array camera (IRAQ |Fazio et al. 2004), die multi- t ure" in mid-infrared spectra of Herbig Ae/Be and T Tauri stars 

band imaging photometer for Spitzer (MIPS, Rieke et al. 2004j), dCohen & WitteboTr] fl98l iNatta et all l2000t IBouwman et al.1 

and the infrared spectrograph (IRS, | HoucketalJ |200J). The | 2 Q01[ l2003b Ivan Boekel et al.ll2004 120051) . This feature is emit- 

Spitzer data allow characterization of the disks around low- ted by si]icate grains with sizes of up t0 several microns 

mass stars up to radii of -10 AU. The observed disk frequen- th at reside in the optical l y thin surface layers of th e disk 

dMen'shchikov & Henningl fl997b IChiang & Goldreichl fl997h . 

* Based on observations performed at ESO's La Silla-Paranal obser- The shape and the strength of the 10/im silicate feature 

vatory under programme 076.C-0470 were found to be correlated and this was interpreted as 
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evidence for grain growth in disks ( van Boek el et al 1 12001 
iKessler-Silacci et al.ll2006t iBouwman et al.ll2008l) . The Spitzer 
IRS has provided spectra of many young objects, covering the 
wavelength range from 5.3 to 38/mi. Analysis of the spectral 
features in the Spitzer IRS spectra provides co nstraints on the 
chemical composition and grain sizes (see, e.g.. lHenning||20l61 
for a recent review). The dust properties derived from Spitzer 
IRS spectra of the members of the MBM 12 and 77 Cha asso- 
ciations suggest that dust processing in the disks occurs very 
early and that radial mixi ng is not efficient ( Mee us et al . 2009; 
ISicilia-Aguilar et al.ll2009l) . 

Because the 77 Cha association is much older than MBM 12, 
it is highly desirable to study an association at an intermedi- 
ate age in order to draw a more complete picture of disk evo- 
lution in sparse stellar associ ations. With an estimated age of 
3-5 Myr (Feig elson et alj 120031) . the eCha association is well 
suited for this purpose. Located at a distance of ~114p c, the 
eCha association was discovered by Feigels on et al .(2003), who 
identified nine members. Th e most massive member, eCha AB, 
has a spectral type of B9. Luhman (2004) found three new 
members of eCha from a survey covering a region with a ra- 
dius of 0°5 around eChaAB. At 3-5 Myr, the e Cha associa- 
tion i s younger than the 77 Cha g r oup (5- 10 Myr, Mamaiek et al. 
I1999L iLuhman & Steeghsll2004l) . lLawson et al.l (120091) provide 
a model-independent way to rank the ages of the nearest PMS 
associations by employing gravity-sensitive spectral features in 
optical spectra. They confirmed that the eCha association is 
younger than the 77 Cha and TW Hya associations. 

In this paper, we will investigate the disks surrounding the 
e Cha members in terms of overall (geometry) evolution and dust 
mineralogy. We arrange this paper as follows. In Section [2] we 
describe the observations and data reduction. In Section [3] we 
present our results, which are then discussed. We summarize our 
efforts in Section |4] 



2. Observations, data reduction, and analysis 

2.1. Targets 

There a re 12 stars reported as the members of eCha in the lit- 
erature (Feigelson et al. l l200llLuhma"nll200l . We collected the 
photometric and spectroscopic data of these stars that are avail- 
able in the literature and list these in Tables [T] [2] and [3] For 
the ten stars whose spectral types have been estimated, we per- 
formed spectral energy distribution (SED) fits in the optical and 
near-infrared spectral range to derive the luminosity of each star 
and its line -of-sight extinction, using the method described in 
Fang et al.l (120091) . We take a model atmosphere spectrum with a 
fixed effective temperature corresponding to the observed spec- 
tral type and generate synthetic photometry with two free pa- 
rameters: the visual extinction Ay and the stellar angular diam- 
eter 9. Both fit parameters are fine tuned to minimize the resid- 
uals between synthetic and observed optical and near-infrared 
photo metry. We adopt a standard extinction law (Cardelli et al. 
119891) with a total-to-selective extinction ratio of Ry = 3.1. We 
used all the available optical photometry and the near-infrared 
J- and H-band photometry. The bolometric luminosity of each 
star is then obtained by integrating the de-reddened model spec- 
trum and adopting a distance of 1 14 pc. For the eCha AB, which 
is a double system, we estimate the luminos ity of the A compo - 
nent from its V-band magnitude (5.34 mag. iTorres et al.l l2008). 
which was corrected for binarity. The resulting extinctions and 
luminosities are listed in Table [3] 



2.2. VISIR imaging 

In the eCha association, four sources (HD104237B-E) are lo- 
cated too close to the bright star HD 104237 A to be properly re- 
solved with the Spitzer telescope. To characterize the disk prop- 
erties of these sources, we observed them with the VLT imager 
and spectrometer for mid-infrared (VISIR) on the night of 2006 
March 17 through the SiC filter, which transmits radiation be- 
tween ~10.7 and ~12.9/jm. We employed standard chopping 
and nodding techniques to eliminate high instrumental and at- 
mospheric background emission, with a chop throw of 8"in the 
north-south direction and a nod throw of 8"in the east-west di- 
rection. We assessed the quality of the individual chop half- 
cycles (here referred to as "frames") by fitting a 2D Gaussian 
function to the A component, which is detected at high S/N in 
each frame. We find that the quality of most frames is good, 
meaning the image of the A component is sharp and round. The 
median full width at half maximum of the 75% best frames is 
0733, and we shift-and-add these into our final image with sub- 
pixel accuracy, using shifts determined by cross correlation. 
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Fig.l. VLT/VISIR image of the HD 104237 region taken 
through the SIC filter, on a square root scaling from -0.1% to 
+6.5% of the peak flux of the HD 104237A component. The po- 
sition of the HD 104237D component is indicated with a circle, 
but there is no significant detection of this source. Coordinates 
are offsets in arcseconds with respect to HD 104237 A. 



In Fig. Q] we show our final VISIR image of the HD 104237 
system. The A and E components, corresponding to sources 3 
and 6 in Grad y et ail J2004). respectively, are clearly detected. 
No sources are seen at the po sitions of the C an d D components 
(source numbers 1 and 5 in iGradv et al.ll2004l) . We did not at- 
tempt an absolute flux calibration of our VISIR data, but instead 
performed point-spread function photometry on the E compo- 
nent, using the A component as the reference star. We find an 
A/E flux ratio of 23.4+0.4. To set the absolute flux scale, we 
integrated the Spitzer spectrum of the A component over the 
transmission curve of the SIC filter. Assuming that the spec- 
trum of the E component is iso-photonic, one obtains a flux of 
551 + 10mJy at 12/jm for this source. At the position of the D 
component, we do not detect a source with confidence and de- 
rive a 3<t upper limit of 18mJy at 12/jm, again assuming an 
iso-photonic spectrum. 
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Table 1. Optical and near-infrared photometr y for memb ers in the e Cha association. Column 5: the B-band photometry of member 
ID#2,5 comes from the Tycho-2 Catalogue dFfog et al] 120 00): the fi-band photometry for member ID#10,1 1,12 is obtained with 
their V-band photometry an d B - V colors dLvo et al.ll2008l) . Column 6: the V-band photometry of member I D#1, 6, 7, 8, 9 c omes 
from Feigelso n et alJ (120031) : the V-band photometry of member ID#2,5 comes from the Tycho-2 Cat alogue dFfog et alJ feOOO): the 
V-band photometry of member ID#10,1 1,12 is obtained from its /-band photometr y and V — I colors dLvo et alj|2008l) . Column 7: 
the R-band photometry of member ID#1,6,7,8,9 comes from Feig elson et al.1 ( 20031). where as the /?-band photometry of member 
ID#10,1 1,12 is obtained with its /-ba nd photometry and R — I colors from (Lvo et al1 l2008l) . Column 8: the /-band photometry 
for member ID #1,6,7,8,9 comes fromlFeigelson et al.l (12003b . whereas the /-band photometry for member ID# 10,1 1,12 is fro m the 
DENIS survey (lEpchtein et alJI 1997b . Columns 9, 10, 1 1: the p hotometry for membe r ID#3,4,6,7 is from lGradv et all (l2004b . The 
photometry for other members comes from the 2MASS survey (Skrutskie et al. 2006|). 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 






RA 


DEC 


B 


V 


R 


/ 


J 


H 


K, 


ID 


Name 


(J2000) 


(J2000) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


1 


CXOUJ1 15908.2-781232 


11:59:07.98 


-78:12:32.2 




16.99 


15.57 


13.83 


12.01 


11.45 


11.17 


2 


eChaAB 


11:59:37.53 


-78:13:18.9 


4.86 


4.91 






5.02 


5.04 


4.98 


3 


HD 104237C 


12:00:03.89 


-78:11:31.0 










15.28 


14.85 


14.48 


4 


HD 104237B 


12:00:04.76 


-78:11:34.8 










11.43 


10.27 


9.52 


5 


HD 104237A 


12:00:05.21 


-78:11:34.4 


6.86 


6.62 






5.81 


5.25 


4.59 


6 


HD 104237D 


12:00:08.39 


-78:11:39.2 




14.28 


13.09 


11.62 


10.53 


9.73 


9.67 


7 


HD 104237E 


12:00:09.43 


-78:11:42.2 




12.08 


11.25 


10.28 


9.10 


8.05 


7.70 


g 


USNO-B120144.7-781926 


12:01:44.42 


-78:19:26.8 




17.18 


15.61 


13.72 


11.68 


11.12 


10.78 


9 


CXOU J120152.8-781840 


12:01:52.52 


-78:18:41.4 




16.78 


15.29 


13.52 


11.63 


11.04 


10.77 


II) 


2MASS J12005517-7820296 


12:00:55.17 


-78:20:29.7 


19.61 


17.85 


16.08 


14.00 


11.96 


11.40 


11.01 


1 1 


2MASS J12014343-7835472 


12:01:43.43 


-78:35:47.2 


20.17 


18.55 


17.15 


15.96 


14.36 


13.38 


12.81 


12 


2MASS J12074597-7816064 


12:07:45.98 


-78:16:06.5 


17.68 


16.08 


14.74 


13.11 


11.55 


10.98 


10.67 



Table 2. Infrared photometry for memb ers in the e Cha associ ation. Columns 4,5,6: the photometry comes from Gra dy et al.l (|2004). 
Column 7: the photometry comes from Luhma net al.l (l2010|). Columns 8,9,10 , 11,12,13: the photomet ry is from the AKARI IRC 
Point Source Catalogue and FIS Bright Source Catalogue (llshihara et al. 20101: Yamamuraet alJl20l6T) . Columns 14,15,16,17: the 
photometry is from WISE Preliminary Release Source Catalog (Wri ght et al.112.0 1 oTT ! 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


(17) 




RA 


DEC 


L' 


M 


11.9/mi 


[24] 


S9W 


L18W 


N60 


WIDE-S 


WIDE-L 


N160 


Wl 


W2 


W3 


W4 


ID 


(J2000) 


(J2000) 


(mjy) 


(mly) 


(mJy) 


(mJy) 


(mly) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mag) 


(mag) 


(mag) 


(mag) 


1 


11:59:07.98 


-78:12:32.2 






















11.00 


10.77 


10.58 


>8.73 


2 


11:59:37.53 


-78:13:18.9 








73.68 


561.4 












5.81 


5.12 


5.10 


4.95 


5 


12:00:05.12 


-78:11:34.7 


9795. 


9091. 


14167. 




13479.7 


13000.1 


10036.2 


8600.1 


5116.1 


4531.8 


4.01 


3.03 


0.81 


-0.71 


6 


12:00:08.30 


-78:11:39.6 


45. 




























7 
8 


12:00:09.32 
12:01:44.42 


-78:11:42.5 
-78:19:26.8 


481. 


329. 


498. 


1055.2 
15.39 














10.36 


9.83 


8.34 


6.74 


9 


12:01:52.52 


-78:18:41.4 








0.80 














10.60 


10.35 


10.04 


>8.67 


10 


12:00:55.17 


-78:20:29.7 








15.54 














10.64 


10.17 


8.48 


6.62 


11 


12:01:43.43 


-78:35:47.2 








55.38 














12.37 


11.59 


7.55 


5.26 


12 


12:07:45.98 


-78:16:06.5 








0.88 














10.51 


10.34 


10.12 


>8.72 



2.3. Spitzer IRS spectroscopy 

The eCha members were observed with the Spitzer Space 
Telescope as part of a large program aimed at studying the 
evolution of circumstellar disks in nearby associations (GO 
proposal 20691, PI: Bouwman). We obtained 7-35 /mi low- 
resolution (R=60-120) spectra with the IRS. The extracted spec- 
tra are based on the droopres products processed through the 
S18.7.0 version of the Spitzer dat a pipeline. Partially based on 
the SMART software package (Higd on et ail 12004). our data 
were further processed using spectral extraction tools developed 
by the Formation and Evolution of Pla netary Systems (FEPS ) 
Spitzer science legacy team (see also iBouwman et al.l 12008). 
The spectra were extracted using a 6.0 pixel and 5.0 pixel fixed- 
width aperture in the spatial dimension for observations with the 
first order of the short- (7.5-14/mi) and the long-wavelength 
(14-35 /mi) modules, respectively. The background was sub- 
tracted using associated pairs of imaged spectra from the two 
nod positions along the slit, also eliminating stray light contam- 
ination and anomalous dark currents. Pixels tagged by the data 
pipeline as being "bad" were replaced with a value interpolated 
from an 8-pixel perimeter surrounding the errant pixel. Low- 
level fringing at wavelengths >20/mi w as removed using the 
irsfinge package (lLahuis & Boogerd2003h . To remove any effect 



of pointing offsets, we matched orders based on the point-spread 
function of the IRS instrument, thereby correcting for possible 
flux losses. The spectra are calibrated using a spectral response 
function derived from multiple IRS spectra of the calibration star 
77 1 Doradus and a MARCS stellar model provided by the Spitzer 
Science Center. The spectra of the calibration target were ex- 
tracted in the same way as our science targets. The relative er- 
rors between spectral channels within one order are dominated 
by the noise in each channel and not by the calibration. We es- 
timate the relative flux calibration within a spectral order to be 
good to wl% and the absolute calibration between different or- 
ders to be accurate to «3%, the uncertainty being dominated by 
uncertainties in the scaling of the MARCS model. 



We use the two-layer temperature distri bution (TLTD ) spec- 
tral d ecomposition routines developed by lJuhasz et al.l ([2009, 
2010) to derive the mineralogical composition of the dust from 
the observed IRS spectra. The TLTD method is applicable to 
disks at most inclinations but do es not work wel l for objects 
with a near edge-on inclination (lJuhasz et al.ll2009b . Hence, we 
can apply it to all sources in this survey, with the exception of 
2MASS J12014343 (see Sect. [3T23T >. The spectrum emitted by 
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Table 3. Stellar parameters for members in the eCha association. Columns 4,5: the values come from iLuhmanl (120041). Columns 
8,9: the masses and ages of the sources ID#1,4,6,7,8,9,10,1 1,12 are estimated from iDotter et al.1 (120081) and Bara ffe et all (fl998). 
respectively. For the sources ID#2A,5 with masses larger than 1 .4 M Q , their masse s and ages are estima ted from Pott er et all (|2008), 
and for the source ID#3 with a mass less than 0.1 M , i ts mass and age are fro m lBaraffe et al.l dl998l) . Column 10: The equivalent 
widths (EW) of Hot line for ID#1, 2, 5,6,7,8,9 come fr om lFeigelson et al.l (l2003l). The EWs of Ha line for ID#10,1 1,12 are estimated 
from the optial spectra presented in lLvo et"ai] d2008). Negative values mean Ha line in emission. Column 1 1: W for weak T-Tauri 
star (WTTS), C for classical T-Tauri star (CTTS), and H for Herbig Ae/Be star. 



(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 




RA 


DEC 




T^g 




Ay 


mass 


age 


Ha 










fs 




ID 


(J2000) 


(J2000) 


Spt 


(K) 


(LO) 


(mag) 


(M ) 


(Myr) 


(A) 


Class 


(M G yr') 


Disk 


(mas/yr) 


(mas/yr) 


Member 


1 


11:59:07.98 


-78:12:32.2 


M4.75 


3161 


0.027 





0.17/0.14 


9.6/7.1 


-6.2 


W 




N 


-35.7±13.6 


-5.8±13.6 


Y 


2A 


11:59:37.53 


-78:13:18.9 


B9 


10500 


99.7 





2.87/... 


2.7/... 


+13 






N 


-35.1±1.5 


4.0±1.9 


Y 


3 


12:00:03.60 


-78:11:31.0 










.../0.015 


.../5 














Y 


4 


12:00:04.00 


-78:11:37.0 








4.5 


0.42/0.17 


4 








Y 






Y 


5 


12:00:05.12 


-78:11:34.7 


A7.75 


7648 


42.8 


0.91 


2.53/... 


3.2/... 


-20 


H 


4.17E-07 


Y 


-36.9±1.4 


-5.6±1.8 


Y 


6 


12:00:08.30 


-78:11:39.6 


M3.5 


3342 


0.12 





0.29/0.30 


3.0/3.8 


-3.9 


W 




N 






Y 


7 


12:00:09.32 


-78:11:42.5 


K5.5 


4278 


0.92 


1.0 


0.90/1.33 


2.1/5.9 


-4.5 


c 


8.59E-10 


Y 






Y 


8 


12:01:44.42 


-78:19:26.8 


M5 


3125 


0.036 





0.16/0.14 


5.1/4.5 


-23 


c 


1.49E-11 


Y 


-31.5±7.7 


0.5±7.7 


Y 


9 


12:01:52.52 


-78:18:41.4 


M4.75 


3161 


0.039 





0.18/0.15 


5.9/5.0 


-7.8 


w 




N 


-28±7.7 


-34±7.7 


7 


10 


12:00:55.17 


-78:20:29.7 


M5.75 


3024 


0.029 





0.11/0.10 


3.1/3.5 


-23.7 


c 


1.05E-11 


Y 


-33.6±7.7 


-2.2±7.7 


Y 


11 


12:01:43.43 


-78:35:47.2 


M2.25 


3524 


0.0038 





0.27/0.49 


4 


-192.4 


c 




Y 


-43.1±8.8 


-6.5±8.8 


Y 


12 


12:07:45.98 


-78:16:06.5 


M3.75 


3306 


0.042 





0.27/0.25 


12.5/11.4 


-7.9 


w 




N 


-68.7±7.8 


-17.9±7.7 


7 



the optically thin surface layers of the disk is approximated by 
the expression 

n m r am , n 2n 

F v = F v , cont + 2_ i 2jA,^i,j J -^B v (T)T—dT, (1) 

where the different dust species in the model are denoted with 
subscripts z; for each species, we include different grain sizes, in- 
dicated with j subscripts. K;j denotes the mass absorption coeffi- 
cient of dust species i with grain size j. The abundances of each 
dust component are indicated by Dj j, while B V (T) denotes the 
Planck function and qa represents the exponent of the adopted 
power-law temperature distribution. The disk atmosphere has a 
range of temperatures between the integration boundaries T a , 
and d denotes the distance between the star and Earth. F ViCont 
denotes the continuum flux from the optically thick disk interior 
and is approximated by 

F v , mat = D -±B v (T*)+D\ J —B v (T)T~dT 

—B V (T)T~ dT. (2) 

Here R+ and are the radius and effective temperature of the 
star. The optically thick inner rim and midplane of the disk are 
assumed to have power-law temperature distributions between 
the T T and T m boundaries, with exponents denoted by qr and 
qm. The Do, D\, and D2 parameters denote scaling factors for 
the emission from the stellar photosphere, the hot inner rim, and 
the disk midplane, respectively. 

Our dust model includes amorphous silicates with olivine 
and pyroxene stoichiometry, crystalline silicates forsterite, en- 
statite, and amorphous silica. For each amorphous species, we 
include three grain sizes with radii of 0.1, 1.5, and 6.0/im, and 
for each crystalline species we include only two grain sizes 
with radii of 0. 1 and 1 .5 /mi, since there was no evidence for 
large crystals (6.0//m) in our data. The opacity curves for the 
amorphous species were calculated using MIE theory, and those 
for the crystalline species were derived using a distributio n 
of hollow spheres (DHS) approximation (Juhasz et al. 2009). 
The optical co nstants were taken fromfServoin &Pirioul?T973h 
for forsterite, iDorschner et al.l (Il995h for amorphous silicate 



with o livine and pyr oxene stoichiomet ry, iHenning & Mutschkd 
(1 19971) for silica, and Jag er et ail (Il998l) for enstatite. To estimate 
the uncertainties in the derived abundances, the TLTD routines 
apply a simple and well-established Monte Carlo technique. The 
TLTD routines add normally distributed noise to the observed 
spectra with an amplitude given by the flux uncertainty in each 
spectral channel, thereby generating 100 versions of each spec- 
trum that are all consistent with the data. For each of these spec- 
tra, the TLTD routines then perform the same compositional fit 
as described above, yielding 100 values for each fit parameter. 
The mean is then adopted as the best-fit value for each parame- 
ter, and the standard deviation in the positive and negative direc- 
tions yields the 1 cr uncertainty on each parameter. The fit results 
will be presented and discussed in Sect. 13.2.51 

3. Results and discussion 

3. 1 . Stellar properties 

In this section, we will first confirm the memberships of the 
eCha stars using proper motion data and then estimate the 
masses and ages of the members. 

3.1 .1 . Proper motions 

Proper motions provide kinematic information of stars, which 
is often a very useful discriminant to separate members of a 
clust er or stellar associ ation from unrelated field ob jects (see, 
e.g jRochau et al.ll2010l) . Recently. iRoser et al.l(l2010l) presented 
a new catalog of proper motions based on the International 
Celestial Reference System (ICRS) using a combination of 
USNO-B1.0 and 2MASS astrometry. This catalog provides 
proper motions of eight eCha members, which are shown in 
Fig. [2] and listed in Table [3] It is clear that most stars have very 
similar proper motions, with the exception of the stars ID#9 
and 12, which move in different directions. We calculate the 
weighted mean of the proper motions of the group members, ex- 
cluding the latter two sources, and find an average proper motion 
of /i Q ,=-36.0mas/yr and /ug=- 1.2 mas/yr for the eCha associa- 
tion, corresponding to a projected linear velocity of ~19.5km/s. 
The 1 and 2cr error ellipses are also shown in Fig. [2] Since the 
stars with ID# 9 and 12 have proper motions that both differ 
by >3 cr from the average group value, this casts doubt on their 



4 



M. Fang et al,: Young stars in e Cha and their disks: disk evolution in sparse associations 



20 



cn -10 



i -20 



-30 



-40 



12 




4- 



-60 -40 
/ii a cos((5) (mas/yr) 



-20 



Fig. 2. Proper motions of known stars in th e eCha association 
taken from the catalog by iRoser et al The mean proper 

motion of the members and the 1 and 2 cr ellipses of the stan- 
dard deviation of the distribution are shown with dashed-line and 
dotted-line curves, respectively. 
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Fig. 3. HR diagram of members (filled circles) of eCha. The 
open diamond marks the underluminous object ID#1 1. The grey 
lines are PMS isochrones of 1, 3, 10, an d lOOMyr, and th e 
dark line is the PMS isochrone of 4Myr (iDotter et alj 12008). 
The dashed lines are the PMS evolutionary tracks for individ- 
ual masses. The open circles are for stars ID#9 and 12, which 
have proper motions that deviate from the other stars. 



membership to the eCha association. No proper motion esti- 
mates are available for the sources with ID# 3, 4, 6, and 7. Since 
these sources show X-ray emission and are spatially associated 
with HD 104237A, it is highly likely that they are eCha mem- 
bers. In summary, we confirm ten of the previously identified 
members as members of eCha, and challenge the membership 
of two stars, ID#9 and 12. The refined memberships of eCha 
are also listed in Table [3] We must stress that we cannot ex- 
clude with certainty that sources ID#9 and 12 are members of 
the eCha association. Both sources may be in binary systems, 
where their orbital motions can alter their proper motions, or 



they may have been ejected from binary systems. Moreover, the 
strength of Li absorption lines in these two stars is comparable 
to other eCha stars, as are the Nal and KI strengths, indicating 
that the t wo stars have ages similar to those o f the other eCha 
members (jEuh man 2004; Lawso n et al. 2009). Source ID#9 is 



also an X-ray emitter dFeigelson et al. 2003 ). Thus both sources 



are clearly young stars and may well have formed together with 
the other eCha members. 



3.1.2. The masses and ages of the eCha members 

Figure|3]shows the HR diagram of the sources in the eCha asso- 
ciation, where the temperatures were adopted from the observed 
spectral types dKenvon & Hartmannlll995t iLuhman et al . 2003) 
and the luminosity was determined as described in Sect. I2.ll In 
the figure, the theoretical pre- main sequence (PM S) evolution- 
ary tracks are obtained from IDotter et"aTI (I2008I) . The masses 
and ages were then deduced by comparing the location of each 
object in the HR diagram with the PMS evolutionary tracks 
(Pott er et al.ll2008l) . For the suspicious members ID#9 and 12, 
masses and ages were estimated, assuming that they are at the 
same distance as the eCha association. As a comparison, we 
also estimat ed the masses and ag es using the PMS evolutionary 
tracks from iBaraffe et al.l d!998l) . In general, both PMS evolu- 
tionary tracks give similar masses and ages for most of sources 
in our sample. 
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Fig. 4. Observed photometry (filled circles) and model pho- 
tometry (open circles) for HD 104237C in the J, H, and K s 
bands. Panel (a) shows model photometry for a 13Mj up brown 
dwarf with an age of 3 Myr, panel (b) shows a model of a 
15.2 Mjup object of 5 Myr. The model photometry was adopted 
from Chabrier elatl d2000h . 



The resulting ages and masses for all stars are listed in 
Table [3] The ages of all the confirmed members, with the ex- 
ception of CXOUJ1 15908.2, are ar ound 3-5 Myr, which i s con- 
sistent with the previous estimate (Feigel son et al.l [2003 ). The 
source ID#1 1 appears sub-luminous by a factor ~30 with respect 
to objects of similar spectral type, placing it below the zero-age 
main sequence (ZAMS) and inhibiting a mass and age estimate 
through placement in the HR diagram. We estimated its mass 
from its spectral type (M2.25) by assuming an age of 4 Myr. 
There is no estimate of the spectral type of source HD 104237B 
(ID#4). We estimated its mass from its position in the Jvs.J-H 
color-magnitude diagram, assuming an age of 4 Myr and using 
model J and H magnitudes from the PMS evolutionary tracks 
(Baraffe et al. 1998: lDotter et al.ll2008t) . The results are listed in 
Tabled 

HD104237C, located ~5'.'3 to the northwest of 
HD 104237A, is a sour ce of X-ray emission seen by Chandra 
dFeigelson et al .] 120031) with a near-infrared counterpart 
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Fig. 5. SEDs of apparently diskless stars. The photospheric 
emission level is indicated with a dotted curve in each panel. 
The circles show the photometry in different bands. The filled 
circle show the upper limit derived from our VLT/VISIR imag- 
ing for ID#6. The thick grey solid lines in panels a, b, d and e 
show the Spitzer IRS spectra. 



(iGradv et alJ [2004). At a distance of 114pc, the projected 
distance between HD 104237C and HD 104237A is -604 AU. 
Figure [4] shows the observed ph otometry for HP 1 04237C in 
the J, H, and K s bands from IGradv et all {2004). Since no 
spectral-type estimate is available for HD 104237C, we estimate 
its mass by comparing the observed near-infrared photometry 
to model p hotometry for young brown dwarfs at ages of 3 
and 5Myr (Chab rier et alJ 120001) . The best-fit masses for the 
assumed ages and the assumed distance of 114pc are 13.0 and 
15.2 Jupiter masses (Mj llp ), respectively. We will adopt these 
numbers as the mass range for HD 104237C, putting this object 
at the boundary between a very low-mass brown dwarf and a 
"free-floating planet". Of the currently known members of the 
eCha association, HD 104237C has by far the lowest mass. 

3.2. Disk properties 

In this section, we first estimate the rates at which material is 
accreted onto the central stars in the eCha member systems 
and then characterize their disks in terms of evolutionary state. 
Finally, we derive the dust properties of those protoplanetary 
disks that show emission from silicate dust in their IRS spectra 
and compare them with the dust properties found for protoplan- 
etary disks in other sparse stellar associations or star-formation 
regions. 

3.2.1. Accretion 

We used the luminosity of the Ha emission line as a proxy for 
the accretion rate. The line luminosity was obtained by integrat- 



ing over the line profile, adopting the best-fit model atmosphere 
spectrum (see Sect. I2.ll ) as the continuum level. We estimated 
the accretion rates from the observed Ha emission-line luminos- 
ity using the empirical relation be tween the latter q uantity and 
the accretion luminosity derived in Fan g et alJ d2009t) : 

log(L. dcc /L e ) = (2.27 + 0.23) + (1.25 ± 0.07) x log(L Ha /L Q ). (3) 

The inferred accretion luminosities are then converted into 
mass accretion rates using the following relation: 



GM*(1 



(4) 



where R m denote s the truncation radiu s of the disk, which is 
taken to be 5R+ dGullbring et al.lll998l) . G is the gravitational 
constant, M* is the stellar mass as estimated from the location of 
each star in the HR diagram, and is the stellar radius derived 
using the SED fitting procedure described in Sect. 12.11 There are 
four sources in e Cha showing signs of active accretion, and their 
accretion rates are listed in Table [3] 

3.2.2. Demographics of disk population in eCha 

In Figs. [5] |6l and |7Jb), we show the SEDs of the ten sources 
for which spectral types are available, together with the best- 
fit stellar model atmospheres (see Sect. 12. U . In total, we de- 
tected nine sources with the Spitzer IRS. Among these, seven 
were confirmed as members of the eCha association, based on 
common proper motion in Sect. l3.1.T1 Two of the members, 
CXOU J 115908.2 (ID#1) and eChaAB (ID#2), do not show 
significant excess emission over the stellar photosphere in the 
IRS spectral range. Also, ID#9 and 12, whose memberships 
were challenged based on their proper motions, do not show 
infrared excess emission in the IRS spectra, strongly suggest- 
ing that they do not (any longer) harbor circumstellar disks. 
For the sources that were not detected with the IRS, we char- 
acterized their disks based on infrared photometric data only. 
The source HD 104237D (ID#6) cannot be spatially separated 
from HD 104237 A with the IRS and MIPS instruments due to 
the comparatively poor spatial resolution of the Spitzer Space 
Telescope at long wavelengths. HD 104237D is not detected 
in our VISIR imaging at 12 /mi, and we derived a 3<x upper 
limit of ~18mJy for its flux at this wavelength. This upper 
limit is just above the expected photospheric emission level, 
as extrapolated from shorter wavelengths and shown in Fig. [5] 
indicating that HD 104237D no longer has a disk. However, 
based on currently available data, the possibility of an evolved 
disk with a large inner hole that would be observable only at 
longer wavelengths cannot be excluded. The underluminous ob- 
ject 2MASS J12014343 (ID#11) shows a very flat SED, which 
is similar to the SEDs of similarly unde rluminous object s dis- 
covered in other star-forming regions (see lFang et al.1 12009'). For 
sources HD 104237B (ID#4) and HD 104237C (ID#3), we only 
have photometry in JHK S bands. The J-Hvs. H-K s color-color 
diagram shows only HD 104237B to still have near-infrared ex- 
cess emission, ind icative of a hot, optically thick inner disk 
(IGradv et al.ll2004h . 

In Fig. Sa), we show the SED of HD 104237 A, which ex- 
hibits strong excess emission at wavelengths beyond ~2/mi. 
HD 104237 A therefore shows clear evidence for a hot, opti- 
cally thick inner disk. We use the 2D RADMC code from 
i ' I j i 

Dullemond & Dominik (2004) to calculate two SEDs, one for a 
flaring disk model and the other for a flattened disk model, given 
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Fig. 6. SEDs of the disk population in e Cha. The photospheric emission level is indicated with a grey dotted curve in each panel. 
The circles show the photometry in various bands. The filled circles show the photometry from our VLT/VISIR imaging. The solid 
line shows the IRS spectrum in each panel. The dashed line in each panel represents the calculated SED of a flaring disk model, and 
the dash-dotted line shows the calculated SED of a flattened disk model. 



the effective temperature (r e ff), stellar radius (/?*), and stellar 
mass (Mi,) of HD 104237 A. In both models, we assume a pres- 
sure scale height (Hp) that varies as a power law with the radius 
(R), Hp/R = R l l 7 . The inner disk radius (R m ) is set at the dust 
destruction radius of T^ust ~1500K, and the outer disk radius 
(Roui) is 500 AU. The scale height at the outer disk radius is set 
as H out /R oul =0.2 for the flaring disk model and as H ou JR out =0. 1 
for the settled/flattened disk model. We use a power-law size dis- 
tribution with an exponent -3.5 for dust grain with a minimum 
grain size of 0. 1 pun and a maximum grain size of 1000 fan. Two 
populations of amorphous dust grains (25% carbon and 75% 
silicate) have been included in the models. In the flaring disk 
model, the disk mass is set to be 0.02 M*, and in the flattened 
disk model, the disk mass (Mdisk) is 0.001 M*. The gas-to-dust 
ratio is taken to be 100 in the models. The disk surface density 
(E) is estimated from the total disk mass, assuming 2 oc R 1 . In 
Fig- Ufa), we show the two model SEDs for the disk inclination 
angle ~ 55°. The SED of HD 104237A is more similar to the 
SED for the flattened disk model than to the SED for a flaring 
disk model. 

In Fig. H|b), (c), and (d), we show the SEDs of group 
members HD 104237E (ID#7), USNO-B 120144.7 (ID#8), and 
2MASSJ12005517 (ID#10), as well as two model SEDs: one 
for a flaring disk and one for a flattened disk model for each 
object. In the calculation, we use the same parameters in the 
disk models as those for HD 104237 A, except that the outer 
disk radius is set to be 200 AU. The actual stellar parameters of 
each object are used in the calculation. Overall, the SEDs of the 
three objects HD 104237E (ID#7), USNO-B120144.7 (ID#8), 
and 2MASSJ12005517 (ID#10) are between the two model 
SEDs. It can also be seen that the excess emission in the MIPS 
24 jum band of source ID#7 is similar to the model SED of the 
flaring disk, whereas its infrared excess at shorter wavelengths 
is substantially lower and closer to the SED of a flattened disk. 
This may indicate that source ID#7 has started to dissipate its in- 



ner disk, while its outer disk remains relatively intact. In contrast 
to object ID#7, the very low-mass objects ID#8 and 10 show a 
fairly uniform depletion of the SED. Their SEDs are similar to 
the SEDs of the flattened disk models. 

Spitzer observations of clusters with ages of sev- 
eral Myrs, such as IC348, NGC2362, 77 Cha, and the 
Coronet cluster, suggest that, qualitatively, two evolution- 
ary paths exist for going from a primordial to a debris 
disk configuration (lLada et al.ll2006l; [sjcilia-Aguilar et al. 2008; 
I Currie & Kenvonll2009t ICurrie et all 20091: ISicilia-Aguilar et alj 
120091: IC urrie & Sicilia-A guilar 2011). Each path shows a char- 
acteristic behavior of the SED: (1) some objects show little or 
no excess emission in the shorter IRAC bands and strong ex- 
cess emission at 24 yum, suggesting that disks are dissipated in an 
inside-out fashion; (2) some objects show infrared excess of an 
approximately uniformly reduced magnitude compared to pri- 
mordial disks over all wavelengths out to 24/im, suggesting a 
reduction of the effective disk height, i.e., the height above the 
disk midplane where the disk becomes optically thick to the stel- 
lar radiation. A reduced disk height causes a smaller fraction of 
the stellar energy to be absorbed and reprocessed by the disk and 
thereby yields an infrared excess of reduced magnitude. This 
may occur if dust coagulation takes place in the disk, causing 
the grains to couple somewhat less well to the gas and allowing 
them to settle towards the midplane. How effective this is also 
depends on the level of turbulence in the disk; grains settle more 
easily in a disk with low turbulence. In some cases, the disks 
could have lost a large number of small dust grains, and become 
optically thin. This type of d isk can also show a significant de- 
pletion of infrared emission (ICurrie & Sicil ia-Aguilarl [201 lh . In 
eCha, source ID#7 shows an SED reminiscent of the first type of 
disk evolution (inside-out fashion), whereas sources ID#8 and 10 
show an SED corresponding to the second scenario (globally de- 
pleted). Both latter sources show strong evidence of dust growth, 
which seems to have progressed particularly far in source ID#8 
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judging from the absence of silicate features at ~10 and 20 /mi. 
The absence indicates that there is no substantial population of 
dust grains with sizes of <10/mi left in the atmosphere of this 
disk. It is interesting to note that source ID#7, which shows evi- 
dence of an inside-out disk clearing, has a silicate feature emitted 
by predominantly small, sub-micron sized grains. 

3.2.3. 2MASS J1 201 4343: a disk system with high 
inclination? 

When the HR diagram of the eCha sample (Fig. [3} is inspected, 
one object stands out in the sense that it lies well below the 
ZAMS. Compared to ~4Myr old stars with the same spec- 
tral type, 2MASS J 120 14343 (ID#11) appears underluminous 
by a factor of -30. Since evidence of group membership is 
strong (optical emission lines, infrared excess, and in partic- 
ular a proper motion in agreement with the eCha group, see 
Sect. 13.1. Il l its underluminosity cannot be simply explained due 
to ID#1 1 being a field star. To understand its nature we should 
resort to other explanations. In Fig. 13a), we show the optical 
spectrum of ID#11, on which Balmer emission lines and for- 
bidden oxygen emission lines with very high EWs are seen. 
Source ID#1 1 also shows a particularly strong infrared excess 
(Fig. 0b))- Objects with these characteristics have been found 
in other star-formin g regions, e.g., Lupus 3 dark cloud, Taurus, 
L1630. and L1641 (IComeron et al.l)200l IWhite & Hillenbrandl 
I2004t IFang et al.ll2009t iMerfn et al.ll2Q10t) . One possible expla- 
nation for this phenomenon is that these sources are systems har- 
boring flared disks with moderately high inclinations where the 
stellar photospheric light is largely absorbed by the material in 
the cold, flared outer disk. We still receive photospheric light, 
but a large fraction of the light we see is scattered off the disk 
surface and has a much reduced total flux. The optical emis- 
sion lines of large EWs may arise in an outflow or disk wind. 
They need not intrinsically be brighter than in similar objects 
with "normal" apparent luminosities; it is the reduced strength 
of the continuum flux, not the intrinsic line strength, that causes 
the EWs to be high. This scenario only works if the line-forming 
region is much larger than the central star, so that at least part 
of the line flux reaches us relatively unhindered, while the pho- 
tospheric continuum is strongly absorbed. In particular, this sce- 
nario also explains why some of these apparently underluminous 
objects show some emission lines like HeI5876A and the Call 
near-infrared triplet (84 98, 8542, and 8662 A) with EWs that ap - 
pear not enhanced (see IComeron et aO2003t IFang et al.l l2009): 
these lines are mainly f ormed in the magnetospheric infall flows 
dMuzerolle et alj [1998 ), which are close to the stellar surface 
and shou l d be a s much occulted as the photospheric continuum. 
iLuhmanl (120041) suggested that 2MASS J12014343 harbors an 
edge-on disk and seen in scattered light. The 10 /mi silicate fea- 
ture of 2MASS J12014343 is seen in emission, suggesting that 
its inclination is substantially different from 90 degrees (edge 
on). 

In Fig. [7] we also compare the optical spectrum and SED 
of source ID#11 with an object that shows very similar char- 
acteristics and was recently found in the Lynds 1 641 cloud 
in Or ion (ID#122, named L1641#122 hereafter, see Fan g et at] 
2009). L1641#122 is also apparently underluminous compared 
to stars of similar spectral type in the same cloud by a factor 
of ~30, similar to ID#11. But contrary to ID#11, it shows the 
HeI5876A, [SII], and Call near-infrared triplet (8498, 8542, 
and 8662 A) in emission, indicating active accretion. The optical 
spectrum of ID#1 1 did not show signs of active accretion at the 
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Fig. 8. SED of the source ID#11. The open circles show the 
broad-band photometry of ID#1 1, and the IRS spectrum of this 
source is shown in black. The dashed lines marked by numbers 
1-5 show model SEDs for a same model disk with different incli- 
nation angles. From model 1 to model 5, the inclination angles 
are 83°, 84°, 84.5°, 85°, and 86°, respectively. In the model, we 
added a few small dust grains to the hole of the disk. The model 
with disk inclination -84.5° best fits the SED of ID#11. The 
dash-dotted line shows the model SED for a disk (inclination 
angle~84°) without small dust grains in the disk hole. The model 
SED generally fits the shape of the observed SED of ID#1 1, but 
fails to reproduce the observed silicate features. 



time of observation. In Fig. |7jb), we compare the median SED 
of the distribut ed disk populatio n in LI 641 (LI 64 ID) and that 
of L1641#122 (IFang et alJl2009h with the SED of ID#1 1. From 
this comparison we can see that the infrared SEDs of ID#1 1 and 
L1641#122 are remarkably similar and that both underluminous 
objects show much stronger mid-infrared excesses than the av- 
erage source in L1641, even though the L1641D population is 
much younger (~1 Myr) than eCha. At wavelengths shorter than 
~8 /mi, the SED of ID#1 1 shows a rapidly decreasing trend. This 
is different from the SED of L1641#122, which shows relatively 
flat SED till near-infrared wavelengths. It suggests that the inner 
disk around ID#1 1 may have been dissipated, consistent with its 
very low or no accretion activities. 

As discussed above, based on the SED of ID#11, the cir- 
cumstellar disk of ID#1 1 may have two properties: (1) the incli- 
nation angle of the disk should be high but less than edge on, 
(2) the inner region of the disk may be evolving. To explore the 
possibility that such types of disks can produce SEDs similar to 
the one of ID#11, we per formed a detailed SED model ing us- 
ing the 2D RADMC code (iDullemond & Dominikjl2004l) . In the 
calculation, the stellar radius and mass for ID#1 1 is given as a 
4 Myr PMS star w ith a spectral type o f M2.25, which is ~1 R e 
and -0.3 M from iDotter et al.1 d2008). By varying M disk , R ouU 
Ri a , and disk inclination angle, we found that the shape of the 
SED of ID#1 1 can be generally fitted by the SED of a disk with 
R in -0.5 AU, M disk -0.0001 AT*, the inclination angle -84°, and 
Hout/Rout ~0.1 (see Fig. [8j. The dust populations in the disk 
are same as those in the disk models described in Sect. 13.2.21 
However, the SED of such a disk model fails to reproduce the 
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Fig. 7. Optical spectrum and SED of 2MASS J12014343 (ID#1 1). (a ) The optical spe ctrum of source ID#1 1 (black line) and of a 
similarly underluminous object in L1641 (L1641-ID#122, grey line, iFang et al.l l2009). (b) The SED of source ID#11. The photo- 
spheric emission level is indicated with a grey solid curve. The open circles show the broad-band photometry of ID#11, and the 
IRS spectrum of this source is shown in black. The photometry of underluminous source L1641-ID#122 is depicted with the open 
squares connected by dashed line for comparison and shows striking resemblance to tha t of ID#1 1 at m id-infrared wavelengths. The 
dash-dotted line presents the median SED of the distributed disk population in L1641 dFang et al.ll2009h . 



observed 10 and 20 ^m silicate features for ID#11. In order to 
produce these silicate features, we added a few small dust grains 
in the disk hole. The dust grains are located close to the inner 
rim of the disk. The grain size distribution is a power-law dis- 
tribution with an exponent -3.5 for dust grain with a minimum 
grain size 0.005 //m and a maximum grain size 0.5 /im. We fixed 
the M disk -0.0001 M+ and H ont /R oul -0.1. By varying the R in , 
inclination angle, the mass of small dust grains, and the location 
of these small dust grains, we found that the SED of ID#1 1 can 
be fitted when R{ n ~0.7 AU, the inclination angle -84.5°, and 
the total mass of the small grains - 10~ n M , which are located 
between 0.6 and 0.7 AU from the central star. 

3.2.4. The disk and accretor frequencies in eCha 

In order to estimate the disk frequency in eCha, i.e., the frac- 
tion of sources that shows evidence of a circumstellar disk in 
the form of infrared excess emission, one needs to identify all 
members of the e Cha association and characterize their infrared 
SEDs. This is a challenging task because the eCha association is 
spread over a large area of the sky and a full inventory is beyond 
the scope of the current work. Instead, we restrict ourselves to 
the region within a ra dius of Of 5 of e C ha AB. This part of the sky 
has been surveyed by Luhma^ d2004l) . and our knowledge of the 
association can be considered essentially complete. In Figs. |5j 
|6l and |7fb), we show the SEDs of ten members with published 
spectral types. Among these, five objects show excess emission 
at near- or mid-infrared wavelengths and thus are harboring cir- 
cumstellar disks. We also include HD 104237 B, which show s 
excess emission at Ks bands according to iGradv et"al] ((2004), 
and obtain a disk frequency among the 1 1 stellar members of 
eCha association of 55*}?% (6/1 lfl This is quite high for a 3— 
5 Myr old population, though this value is clearly subject to low- 
number statistics and it is premature to conclude that the disk 
frequency in eCha is high based on these data alone. However, 
as we will see later in this section, the five other sparse associ- 
ations of which we estimated the disk frequencies, Taurus, the 
Coronet cluster (CrA), MBM 12, 77 Cha, and TWHya, show a 
similar trend. The ensemble of data provides strong evidence 

1 The uncertainties in disk freq uencies for small samp les are all esti- 
mated in the manner described by Buraasser et al.l d2003h . 



that disk lifetimes in sparse associations are longer than those 
in more crowded environments. 

In Fig. [9] we illustrate this by showing the estimated disk 
frequencies in the aforementioned sparse associations and a 
number of other star-formation regions as a function of their 
age (see Appendix lAl for a detailed description). The disk fre- 
quencies are estimated using only IRAC data or infrared data 
at wavelen gths shorter than ~8 ^m, and the ages are all esti- 
mated from Bara ffeetal.l (fl998). Figure [9] also shows a fit to the 
observed disk frequencies of all regions, except the six sparse 
associations, of the form /disk=exp(-f/To). We find a value of 
to = 2.8+0.1 Myr to yield a good fit (reduced x 1 - 1-8) to the over- 
all distribjrdon. This a grees with the value of -3.0Myr found 
by iFedele et al.1 (I2010J) in an earlier, similar study. Because the 
disk lifetime in any give n environment may depend on the mass 
of the central star (e.g.. iKennedv & Kenvonll2009l) it is useful 
to investigate the low-mass population separately. In Fig. |9|b), 
we show the disk frequency among stars with estimated masses 
in the 0.1 -0.6 M Q range for the same star-forming regions as 
shown in Fig. [9] but including only those for which the low- 
mass population has been characterized. We find essentially 
the same typical disk dispersal time as for the whole mass 
range (to = 2.5+0.2 Myr) (reduced ^ 2 = 1.1). Again, the sparse 
associations are deviant in the sense that they show systemat- 
ically higher disk frequencies. We fitted the disk frequencies 
in the sparse associations of the form /disk=exp(-?/To), which 
gives to = 4.3+0.3Myr (reduced x 2 -^- 6) among all sources and 
to = 4. 1+0.5 Myr (reduced^ 2 =1.5) among sources with masses 
in the 0.1 -0.6 M range. The estimated lifetimes of disks in the 
sparse associations are longer than those in the dense (cluster) 
environments with more than 2cr confidence. As noted in Fig. [9] 
disk frequencies for the sparse associations are nearly constant 
at 60-70% for ages below 4 Myr, and decrease slowly at higher 
ages. This is suggestive of abimodal mechanism responsible for 
disk dissipation in these regions. Very early in the evolution of 
the associations, 30^-0% of the disks may have been dissipated 
quickly by some efficient mechanism, fo r example, interaction 
with binary companions in close systems dBouwman et al.l2 006: 
iKraus et aLll2012l) . After the fast disk-dissipation phase, the re- 
maining disks would evolve slowly via other mechanisms, e.g., 
viscous evolution, planet formation, or photoevaporation. 
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To estimate the accretor frequency, i.e., the fraction of mem- 
bers showing signs of active accretion, we turn to the Ha line 
as a diagnostic. Estimates of the Ha EWs are available for eight 
of the eCha members, of which five a r e cla ssified as accretors 
according to the criteria in Fan g et al.l d2009l) . This brings our 
estimate of the accretor frequency in 6 C ha to 63 + |g% , aga in ob- 
viously subject to low-number statistics. Fedele et al. (201(J) in- 
vestigated the age dependency of accretor frequency and derive 
an empirical relation /acc=<?~'/ 2 ' 3 , where t is in units of Myrs, by 
fitting an exponential profile to the observed accretor frequen- 
cies of a number of star-formation regions. For eCha, aged 3— 
5 Myr, the empirical relation predicts an accretor frequency of 
~10-30%, substantially below the observed value. 

In conclusion, we can state that the disks around stars in 
sparse associations evolve more slowly than those in denser en- 
vironments. 




For 0.1-0.6 M Q '. 




2 4 6 8 10 12 14 
Age (Myr) 



Fig. 9. Disk frequencies for different clusters/star-formation re- 
gions j3]2tted_as_afimcJion of their ages, which are all estimated 
from lBaraffeetalJ (I1998I) (see Table EOV In panel (a) we plot 
the disk frequency among all known members, in panel (b) we 
show only the low-mass (0. 1-0.6 M ) population. The filled pen- 
tagram represents the sparse stellar associations, Taurus, CrA, 
MBM 12, eCha, //Cha, and TW Hya in panel (a), and the sparse 
stellar associations, Taurus, MBM 12, eCha, 77 Cha, and TW Hya 
in panel (b). The open circles show the disk frequencies in the 
compact clusters or the OB associations (see Appendix [A] for 
a detailed description). In each panel, the dashed line repre- 
sents the best-fit exponential decay for all regions, excluding the 
sparse associations and the dash-dotted line represents the best- 
fit exponential decay for the sparse associations. 
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Fig. 10. The 10/um silicate features fitted using the TLTD 
method (solid red lines). The observed spectra are represented 
as filled circles with the errors in grey. 
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Fig. 11. IRS spectra in the 17-37 /jm range and spectral fits using 
the TLTD method. The observed spectra are represented as filled 
circles with the errors in grey, whereas the fits are shown in red. 

3.2.5. The dust properties of disks 

We have used the TLTD spe ctral decomposition method devel- 
oped by Juhas z et"aT1 (120091) to analyze our IRS spectra of the 
eCha members. In this section, we will first present the derived 
dust properties of the protoplanetary disks in e Cha. Then we will 
join the eCha sample with a collection of cool T Tauri stars in 
the MBM 12 and 77 Cha associations, as well as with the Coronet 
cluster and the sample in the cores-to-disks (c2d) legacy pro- 
gram, for which IRS spectra have been previously analyzed us- 
ing identical or similar methods. The combined data set is then 
used to do a statistical study on the dust properties of protoplan- 
etary disks around the cool stars. 

(a) The IRS spectra of the e Cha disks 

There are five eCha members for which we have IRS spectra 
that show evidence for a protoplanetary disk. Four of these show 
the well-known silicate feature between 8 and 13/im in emis- 
sion (ID#5,7,10,1 1, see Fig. ITOt. Before using TLTD routines to 
derive the dust composition, we divided the IRS spectra into a 
"short" (7-17//m) and a "long" (17-37/mi) wavelength regime, 
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Table 4. Silicate compositions estimated from IRS spectra with TLTD. The values are percentages in mass fraction. 
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which are analyzed independently. Since the short wavelength 
part of the spectrum is dominated by warm disk regions closer to 
the central star than the cooler regions that dominate long wave- 
length data, we can search for radial gradients in the dust prop- 
erties. 

The derived dust properties are listed in Tables |4] and [5] In 
Table |H we present the outcome of the fits in terms of mass 
abundances of the different species. In Table[5] we give the mass- 
averaged grain sizes of the amorphous and crystalline silicates, 
as well as the fractional contribution of the crystalline species 
to the total dust mass present in the optically thin disk atmo- 
sphere. As shown in Table [5] the crystalline fractions of the ob- 
jects in e Cha are comparable t o those of young stars with similar 
spectral type in th e literature dSiciha-Aguilar et al.ll2008l [2009; 
iMeeus et al.ll2009t Uuhasz et al.ll2010l) . Also, in the warm disk 
regions, the dust grains of our objects show typical larger sizes 
than those in the interstellar medium, suggesting dust growth in 
these regions. Based on simulations in a ID vertical column of a 
protoplanetary disk, IZsom et all (1201 ll) suggest that high values 
of turbulence (a ~0.01) are needed to explain the existence of 
grains with sizes of a few microns in disk atmospheres at ages 
of several Myr. 

We now briefly describe the individual objects: 
HD 104237A (ID#5) Uuhasz et al l d2010h fitted the IRS spec- 
trum of this object using amorphous dust with olivine and py- 
roxene stoichiometries, as well as silica with grain sizes of 0.1, 
2.0, and 5.0 fim, crystalline species forsterite and enstatite with 
grain sizes of 0.1 and 2.0/mi, and polycyclic aromatic hydro- 
carbons (PAHs). They divided the IRS spectrum into two sec- 
tions: 5-17/im and 17-35/im, which are similar to those we 
use. From the short wavelength part, they derived mass-averaged 
grain sizes of the amorphous and crystalline silicates, which are 
4.6/itn and lAfaa, respectively, and a mass fraction of 9.5% in 
crystalline silicates. In the long wavelength part of the spectrum, 
they find 0. 1 yum and 0.4 yum for the mass-averaged grain sizes 
and 7.1% for the crystallinity. Even though we use a slightly 
different set of grain sizes in our fit, ou r results agree well 
with those derived by Juhasz et al. (2010), with the exception 
of the mass-averaged grain size of the amorphous silicates de- 
rived from the long wavelength channel. Here we find a some- 
what higher value. This may be due to the fact that the mass 
absorption coefficients of amorphous silicates with sizes of 0. 1 
and 1.5 jjm are very similar in the 17-37/im spectral range. 
2MASS J12005517 (ID#10) We find that the spectrum of this 
object is best reproduced using an extraordinarily high mass 
fraction of crystalline silicates of ~49%. This may be a con- 
trast effect: because of the very low luminosity of ID#10, its 
10/im silicate feature arises mainly in the inner ~0.5AU of 
the disk. The central re gions of disks can be highly crystalline 
( Ivan Boekel et al. 2004), possibly leading to very high apparent 
crystallinities if only the very central disk regions contribute to 
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Fig. 12. (a) IRS spectrum of source ID#8. The positions of the 
C2H2 and HCN rovibrational bands are marked, (b) IRS spec- 
trum of the source ID#10. The dotted vertical lines mark the 
features of enstatite. (c), (d), and (e): The mass absorption co- 
efficients of enstatite and amorphous silicates (with olivine and 
pyroxene stoichiometry) with grain sizes of 0. 1 yum(solid lines), 
1.5 //m(dashed lines), and 6.0 / um(dash-dotted lines). The dotted 
lines in panels (b) (c) are identical. 



the p art of the spectrum used in the mineralogical analysis (see 
also lApai et alll2005l) . In Fig. |T2j we compare the IRS spec- 
trum of ID#10 in the 7.8-14.3 /im spectral range with the mass 
absorption coefficients of enstatite and amorphous silicates with 
different grain sizes. The spectral signature of enstatite is clearly 
present in the spectrum of ID#10. 

2MASS J12014343 (ID#11) As discussed in Sect. 13331 source 
ID#11 may be harboring a disk with a high inclination. 
Therefore, the IRS spectrum of ID#1 1 can be moderately red- 
dened by the cold outer disk. In order to deredden the spectrum, 
we would need to know the magnitude of the extinction and the 
proper extinction law at the mid-infrared wavelengths, both of 
which are not well known. We therefore applied no reddening 
correction, took the observed spectrum at face value, and fed it to 
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Table 5. Fitting results for IRS sepctra with TLTD. Columns 2, 
3, 4: the average sizes of amorphous ((a)* nsil ) and crystalline 
grains ((a)^, stsil ) and the mass fractions of crystalline grains 
(f* yst ) are derived from 7-17/im. Columns 5, 6, 7: the average 
sizes of amorphous ((a)^ msil ) and crystalline grains ((a)™ sts ji) 
and the mass fractions of crystalline grains (f™ rt ) are derived 
from 17-35 /um 
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the TLTD routines. Since the TLTD routines only work properly 
for low-extinction sources, our results for source ID#1 1 should 
be regarded as tentative. 

USNO-B120144.7 (ID#8) The luminosity and spectral type of 
this source are similar to those of ID#10. Both sources also show 
similar accretion activity (see Table |3]l and similar SEDs (see 
Fig. [(J- However, in contrast to ID#10, the IRS spectrum of ID#8 
does not show a silicate feature around 10/im (see Fig.fTZt. This 
suggests that the silicates in the warm disk atmosphere of ID#8 
have grown to sizes larger than those in the disk of ID#10. In 
the IRS spectrum of ID#8, we clearly detect the C2H2 rovibra- 
tional band at ~13.7 fj.m, but do not find the HCN rovibrational 
band at ~14.0/um (s ee Fig. [T2b . This is consistent with what 
Pasc ucci et alJ d2009h discovered in Spitzer spectra of other cool 
stars, namely, that the C2H2/HCN flux ratios from cool star disks 
are large with a median of ~3.8, which is an order of magnitude 
higher than the median C2H2/HCN flux ratio in spectra of disks 
surrounding sun-like stars (~0.34). 

(b) Dust properties of disks around cool T Tauri stars 

We now combine the dust properties derived for the eCha 
members with those derived for several other regions and do 
a statistical analysis of the combined data set. We collected 
the available literature data for protoplanetary disks surround- 
ing young stars of spectral type M0 or later in the M BM 12 
(Mee us et all l2009t) . 77 Cha dSicilia-Aguilar et al.ll2009h associ- 
ations, the Coronet cluster dSicilia -Aguilar et al. 2008), a nd th e 
sample in the c2d legacy program (iOlofsson et al.l l2009l 12010). 
The dust properties for the targets in MBM 12, 77 Cha, and the 
Coronet cluster are derived from Spitzer IRS spectra using the 
same TLTD spectral decomposition routines that we used for 
the e Cha data. As described in Sect. 13.2.51 the IRS spectra are 
divided into two regions, a short (7-17yum) wavelength part 
tracing the warm inner disk regions and a long (17-37/jm) 
wavelength part that is more sensitive to cooler regions further 
from the central star . For targets in the c2d legacy program, 
IOlofsson et alJ d2010l) performed the spectral decomposition us- 
ing the B2C method, which includes continuum emission and 
two main components, i.e., a warm and a cold one responsi- 
ble for the 10/i/m and 20-30/jm emission features, respectively. 
This more or less corresponds to the short and long wavelength 
regime of the IRS spectra in the TLTD spectral decomposition 
approach. 

In Fig. [13] we compare the dust properties in the warmer 
and cooler regions. We compare the average sizes of dust grains 
for both the amorphous and crystalline grain population, along 



with the mass fractions of crystalline silicates. As shown in 
Fig-G2Ia)(b), the average sizes of the amorphous and crystalline 
silicates in the warmer regions of protoplanetary disks are sub- 
stantially larger than those in the cooler disk. This may indi- 
cate that dust growth has been more efficient closer to the cen- 
tral star than in regions at larger distances. However it may also 
mean that the disks are more turbulent in the warm inner regions 
than further out, which allows larger grains to remain visible in 
the disk surface instead of settling to the midplane. A combina- 
tion of both effects is also plausible. Figure fOl c) shows that the 
dust in the warm inner regions of protoplanetary disks generally 
contains a higher fraction of crystalline material than the dust 
in cooler regions. This is consistent with earlier findings (e.g., 
Ivan Boekel et~al]|2004l:lMeeus et alj|2009l) . 

We estimated the accretion rates for all stars from their Ha 
line luminosity using the method described in Sect. 13.2.11 and 
list these in Table IB. II In Fig. [14] we plot the average grain 
sizes of amorphous silicates against accretion rates. The data 
suggest a positive correlation between the average grain sizes 
of the amorphous silicates and the accretion rates if the latter 
is above ~1O~ 9 M yr -1 (see Fig. [T4l . Below this value, the 
grains in the majority of disks are small, independently of the 
accretion rate. In order to see whether there is any correlation 
between both observables, we apply a Kendall r test. If two 
datasets are fully correlated, this test returns a value of r = 1 . 
If they are anti-correlated, we get r = - 1 , and if they are in- 
dependent, we obtain r = 0. The Kendall r test also returns a 
probability p, which is smaller when the correlation is more 
significant. We use the Kendall t test to evaluate the possible 
correlation between the average grain sizes of amorphous sili- 
cates and the accretion rates, which yields r=0.79 and p=0.006 
for M acc > 1O~ 9 M yr -1 , and r=0.16 and p=0.58 for M acc < 
10~ 9 M Q yr . Thus there is a significant correlation between the 
average grain sizes of amorphous silicates and the accretion rates 
when M acc > 10~ 9 M yr -1 . A possible explanation for this rela- 
tion is that both accretion and the presence of large grains in the 
disk surface require some level of turbulence in the disk. Thus, 
both a large average grain size and a high accretion rate are tale 
tell signs of a turbulent disk, though there need not be a direct 
causal connection between both observables. 

There are three outliers in Fig. [141 MBM 12-10, G-14, and 
Sz76, which are all classified as WTTSs due to their small Ha 
EWs. Yet, all of them show large average grain sizes in their 
IRS spectra. G-14 shows a quite globally depleted SED, which 
can be d ue to strong overall dust se ttling and/or grain growth in 
the disk dSicilia-Aguilar et alJl2008l). MBM12-10 and Sz 76 can 
be cla ssified as transition disks dMeeus et al . 2009; Wahhai ~et alJ 
[2oToh . It is still unknown which mechanisms are responsible for 
the evolution from normal disks to transition disks. One of the 
proposed mechanisms is dust growth and settling, which is be- 
lieve d to work in some of the tran sition disks in the Cep OB2 re- 
gion ( Sici lia-Aguilar et al.ll201 la|). How ever, in the Chamaeleon 
I star-forming region. lManoi et all d201 ll) found that some of the 
transition disks show 10 fim features from relatively unprocessed 
grains compared with normal disks. This indicates that, at least 
in disk surfaces, the dust population has not undergone substan- 
tial grain growth in these objects. Thus that a transition disk 
geometry does not necessarily go hand in hand with large sil- 
icate grains dominating the mid-infrared spectra. Nevertheless, 
the deviant behavior of sources MBM12-10, G-14, and Sz76 
from the correlation between the average grain size of the amor- 
phous grains and the accretion rate is plausibly due to their spe- 
cial disk structures (with inner hole s) or extreme grain growth 
(see also lSicilia-Aguilar et alJl201 Tah . 
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Fig. 13. (a) Comparison of the mass-averaged sizes of amorphous dust grains as derived from the shorter and longer wavelength 
part of the IRS spectra (indicated with "W" and "C " superscripts for "war m" and "cold" disk regions) of sources in e Cha, MBM 12, 
T] Cha, and the Coronet cluster, and the targets from Olofsson et al. (2010). The circles are for the targets in MBM 12, r\ Cha, and the 
Coronet cluster. The diamonds show the targets from lOlofsson et all d2010l) . The pentagram marks object ID#10 and 1 1 in eCha. 
(b) Similar to (a), but for mass-averaged sizes of crystalline dust grains, (c) Comparison of the mass fraction of crystalline dust 
grains as derived from the shorter and longer wavelength part of the IRS spectra. The symbols are similar to panel (a). 



ISicilia-Aeuilar et ail (120071) have discovered a correlation 
between the average grain sizes of silicates and accretion 
rates that is similar to what we find for accretion rates above 
1O _9 M yr . In addition, they argue that turbulence in disks 
prevents large grains from settling into the disk interior, where 
they are invisible to Spitzer, and also promotes accretion, thereby 
leading to a correlation between both observables. What is new 
in the current study is that the correlation between the accre- 
tion rates and grain sizes breaks down for M acc <10~ 9 M Q yr _1 , 
suggesting that the turbulence required to sustain such accretion 
rates is insufficient to support large grains against settling. We 
fit the relation between the average grain size of the amorphous 
dust and the accretion rate with the following curve: 

w _ / 38.2 + 4.2 x LogM acc if M acc > lO" 9 M yr" 1 
Wam.sU. I ^ 1.5 if Macc < 1O" 9 M yr" 1 . 

Here, (a)™ si] (in/im) is the average size of amorphous sili- 
cates in the warmer disk regions, and M acc is the accretion rate 
inM yr -1 . 

In Fig. [14] disks in younger MBM 12 show more processed 
dust grains than those in relatively older r\ Cha. Since younger 
disks tend to show more processed dust grains, we cannot ex- 
plain the positive correlation between (a)™ u as a global effect 
of disk evolution. We would expect older disks to show bigger 
dust grains because dust grains are expected to grow with disk 
evolution. 

Figure Q3] shows the mass fraction of crystalline material in 
the "warmer" and "colder" disk regions versus the effective tem- 
peratures of the central stars for our sample of cool T-Tauri stars. 
We use the Kendall t test to evaluate any possible correlation be- 
tween both quantities, which yields r=-0.25 and p=0.07 for the 
warm disk part (panel a in Fig. [T3J, and r=0.19 and p=0.19 for 
the cooler disk part (panel b in Fig.[T5ll. Thus there is no signifi- 
cant correlation between the mass fraction of crystalline silicates 
in the disk and the stellar effective temperature for cool T Tauri 
stars. 

In Fig. Q2] we plot the crystallinity of the disk material 
against the mass accretion rates. Once again, we run the Kendall 
r test, and find essentially no significant correlation in the warm 
disk regions (r=-0. 27, p=0.10) and the c older parts of the disks 
(t=-0.02, p=0.88). iGlauser et ail d2009t) suggested that irradi- 
ation of dust grains by energetic ions from the stellar winds of 



MBM 12-10 
<-ffi 




-12 -11 -10 -9 -8 

Log 10 M occ (M /yr) 



Fig. 14. The average sizes of amorphous grains plotted versus 
the accretion rate of the central stars. The grain sizes are derived 
in the warmer region of the IRS spectra. The filled symbols show 
the CTTSs, and the open symbols represent the WTTSs. The tar- 
gets in MBM 12 are shown as circles, in 77 Cha as squares, and in 
the Coronet cluster as trian gles. The diamonds show the targets 
from lOlofsson et all d2010l) . The pentagram marks object ID#10 
in £ Cha. The thick dashed lines are the fit to the relation between 
the accretion rates and the grain sizes for M acc >10~ 9 M yr -1 . 

young stars can amorphize the surface layer of the protoplane- 
tary dust very efficiently, thus erasing any correlation between 
crystalline mass fraction and stellar parameters, such as bolo- 
metric luminosity, effective temperature, accretion rate, or disk 
geometry. Abraha m et all (120091) discovered that episodes of in- 
creased accretion may create new crystals. If strongly variable 
accretion is characte ristic of T Tauri stars, as suggested by, e.g. 
Murp hy etail (1201 ll) and Fang et al. (submitted to ApJs), this 
may also account for the lack of obvious correlations between 
any disk or stellar properties and crystallinity of the dust in the 
disk surface. 

4. Summary 

We have studied the members of the e Cha association, focusing 
on the properties of the central stars and their circumstellar disks. 
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Fig. 15. The mass fractions of crystalline material compared to 
the effective temperature of the central stars (which correlates 
closely with the bolometric luminosity, see Fig f3j. The crys- 
tallinities are derived from the 7-17 fim part of the IRS spectra, 
(b) similar to (a), except that the crystallinities are derived from 
in 17-37/mi part of the IRS spectra. The symbols are identical 
to those in Fig. [14] 



tively high disk frequencies. Disk evolutions appears to proceed 
substantially more slowly in sparse associations compared to 
denser environments. In addition, the disk frequencies in sparse 
associations are almost constant at ages less than ~4 Myr. 

The 13.7jum rovibrational band of C2H2 is detected in the 
IRS spectrum of USNO-B 120144.7. 

We derive the mineralogical composition and grain size dis- 
tribution of the (sub-) micron-sized dust in the disk atmosphere 
using the TLTD method. We find that the average grain sizes 
and fractions of crystalline material are higher in the warm in- 
ner disk regions that dominate the short wavelength part of the 
IRS spectra compared to cooler regions at a longer distance from 
the central star that contribute mostly to the longer wavelength 
range covered by the IRS. We also find that the average sizes of 
amorphous grains in the warm inner disk regions show a posi- 
tive correlation with the accretion rates if the latter is higher than 
-lO-'Moyr 1 . 
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Fig. 16. The mass fractions of crystalline silicates plotted ver- 
sus the accretion rates of the central objects. Panel (a) shows the 
crystallinities as derived from the 7-17 /jm part of the IRS spec- 
tra, panel (b) shows the same quantity derived in the 17-37 /xm 
spectral region. The symbols are identical to those in Fig. [14] 
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We used a combination of available archive data, our own Spitzer 
spectroscopy, and VLT/VISIR imaging data. 

Using proper motions, we confirm the membership of 
most studied stars but question the membership of stars 
CXOUJ120152.8 (ID#9) and 2M AS S J 1 2074597 (ID# 1 2) . We 
estimated the masses and ages of the eCha members and find 
HD 104237C to be a sub-stellar object with a very low mass of 
13-15 Mj up , putting it at the boundary between brown dwarfs 
and "free-floating planets". 

The object 2MASSJ1 20 14343 (TD#1 1) is unusually faint at 
optical wavelengths and also shows exceptionally large EWs 
of some optical emission lines and an exceptionally strong in- 
frared excess. Similar objects have been discovered in other 
star-forming regions. We tentatively explain this behavior with 
a flared disk seen at moderately high inclination in which the 
cold outer disk regions cause sufficient extinction to effectively 
screen the central star, which is then seen mostly in scattered 
light, but allow most of the infrared light from the warm disk 
regions to pass. The protoplanetary disks surrounding the cool 
stars USNO-B 120144.7 and 2MASS J12005517 show evidence 
of a reduced height of the optically thick disk due to dust set- 
tling. The disk around HD 104237E shows evidence of partial 
dissipation in its inner part, while its outer disk remains essen- 
tially intact. 

We found that both disk frequency and accretor frequency 
in the e Cha association are higher than those in relatively dense 
clusters of similar age. Five other sparse stellar associations for 
which data are available in the literature also show compara- 



Appendix A: The disk frequencies in star-formation 
regions 

In Table IA.1I we list each star-formation region (SFR) used 
in Fig. [9] In total, 16 SFRs are included with median ages 
ranging from ~1.5 to ~ 13 Myr. Most of the SFRs in the ta- 
bles have an age estimat e using the PMS evolutionary tracks 
from Baraffe end] (1 19981) in the literature. The references for 
these ages are also listed in the table. The 5 SFRs, i.e., e Cha, 
TW Hya, L988e, Tr 37, and NGC 7160, are without ages esti- 
mated from the tracks o f Baraffe et al. (1998). We estimate their 
ages using these tracks with a mixing length parameter l/Hp-l. 
In Table IA.1I we give the disk frequency for each SFR as well 
as the references. We corrected these disk frequencies by remov- 
ing the transitional disk objects that only show excess emission 
at wavelengths longer than 8 //m. For the disk frequencies esti- 
mated in this work, we distinguish the disk po pulation from disk - 
less ones using the SED slope criteria from lLada et al.l (|2006). 
The disk frequencies in the sparse stellar association, MBM 12, 
e Cha, 77 Cha, and TW Hya tend to be systematically above those 
of the other SFRs with similar ages. 

MBM 12 In this association, there are 12 known mem- 
bers with K to M spectral types. Among these, eight ob- 
jects (MBM 12-1,2,3,4,5,6,10,12) were detected by the Spitzer 
IRS, and seven objects show evidence for a circumstellar disk 
(MBM 12-2,3,4,5,6,10,12). Three of t he seven objects a re sug- 
gested to be transition disk candidates dMeeus et al.l2009l) . In the 
recently publishe d catalog from the W ide-field Infrared Survey 
Explorer (WISE, Wrig ht et all 1201 Oh the four objects lacking 
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Fig. A.l. SEDs of two stars in MBM 12 associations. The pho- 
tospheric emission level is indicated with a grey dashed curve in 
each panel. The open circles show the photometry in different 
bands. 
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Fig. A.2. SEDs of two stars in the 77 Cha association. The pho- 
tospheric emission level is indicated with a grey dashed curve in 
each panel. The open circles show the photometry in different 
bands. 



IRS spectra (MBM12-7, 8,9,1 1) are detected in three photo- 
metric bands of WISE: Wl (3.4 /mi), W2 (4.6/mi), and W3 
(12/mi) bands. With the new infrared photometry, we identify 
only MBM12-9 (see Fig. lA. lb as having an IR excess at 12 //m, 
whereas the other three sources have WISE fluxes consistent 
with the stellar photospheric emission. The object MBM 12-8 
shows strong Ha emission (EWhq.=-120A), pointing at active 
accretion. However it shows no excess emission at 12 /mi and 
shorter wavelengths, which indicates that the warm inner disk 
regions have already dissipated. Possibly this is a transition disk 
object with a large opacity hole in the inner disk, where some gas 
is still present. This allows accretion to proceed. Observations at 
longer wavelengths are needed to confirm the transition disk na- 
ture of MBM 12-8. For the time being, we do not consider the 
object to have a confirmed circumstellar disk. Thus, the disk fre- 
quency for all known members in MBM 12 is ~67+}°% (8/12), 
while for members with masses of 0. 1-0.6 M Q , it is ~ 67^7% 
(6/9). 

77 Cha In the 77 Cha association, 18 members have been 
observed with Spit zer. The disk frequency is estimated to 
be_~44+]2% (8/18) (Megeat h et all 12005b ISicilia-Aguilar et all 
2009). The sources RECX-3 and RECX-4 only show excess 
emission at wavelengths longer than 10 /mi. In this work, we 
restrict our disk census to the Spitzer IRAC bands, thus ex- 
cludin g the two sou r ces w hen calculating the disk fraction in 
77 Cha. Murph y et alj ( 120101) proposed that seven candidate low- 
mass members of 77 Cha comprise a halo surrounding the cluster 
core. Two of these may have disks that are actively ac creting, as 
indicated by their Ha emission (Murphy et al. 2010). We have 
found all seven objects in the WISE catalogue. Only one of these 
(J08202975-8003259) shows excess emission in the WISE data 
(see Fig. [A2]i. The source 2MASS J08014860-8058052 does 
not show clear excess emission at wavelengths up to 22 //m (see 
Fig. lA. 21 ). Thus, we presently consider a possible disk around the 
latter object to be confirmed. Including the seven new members, 
the disk frequency of 77 Cha is ~28^'°% (7/25). For the members 
with masses of 0.1-0.6 M Q , the disk frequency is ~32+' 2 (6/19). 



TW Hya To estimate the disk frequency of this associa- 
tion, we first constructed a catalog of members of this group. 
Our catal og is based on th e TW Hya membership criteria as re- 
fined by Mamaiek (2005|), and includes TWA-1, 2A, 2B, 3A, 
3B, 4, 5A, 5B, 6, 7, 8A, 8B, 9A, 9B, 10, 11 A, 11B, 13A, 13B, 
14, 15A, 15B, 16, 20, 21, 23, 25 26, 27, 28, plus new iden ti- 
fied members TWA-29 30A, 30B (lLooper et alil2007l 120 1 OalTbh . 
Among these members, TWA-1, 3A, 4, 7, 1 1A, 27, 28, and 30B 
show evidence fo r disks with excess emission at wavelengths 
shorter than 8 um dJayawardhana et alll999l:lRiaz & Gizisl 2008: 
Loop er et alll2010allbl) . TWA-30A may harbor a disk because its 
optical a nd infrared spectra s how emission lines due to accretion 
activity (lLooper et al.ll2010ah . However, there is no infrared data 
to suggest whether TWA-30A shows excess emission at wave- 
lengths shorter than 8 //m. Thus, we do not count this source 
when calculating the disk fraction in TW Hya. The disk fre- 
quency for all known members is estimated to be ~25^?%(8/32), 
and the disk frequency for members with masses of 0. 1-0.6 M Q 
is ~19+' 3 % (3/16). 

Appendix B: The dust properties around M-type 
PMS stars in other regions 

In Table IbTI we list the targets used in Figs. [13] [14] Q3J and 
[TBI The total luminosity, extinction, and accretion rate listed 
here for each o bject are derived us ing the procedure described in 
Sect. and in Fang et all (|2009). The stellar masses and ages 
are estimat ed by compa r ison t o theoretical PMS evolutionary 
tracks from iDotter et alJ d2008l) . We collected the dust proper- 
ties for each disk from the literatu re (Sicilia-Aguilar et al. 2008; 
i Meeus et all 120091: ISicilia-Aguilar et alj 120091 lOlofsson et all 
12010^ These include the average sizes of amorphous and crys- 
talline silicate grains, along with the mass fractions of crystalline 
silicate grains, and were derived using either the same TLTD 
method that we used for the analysis of the spectra of eCha 
members (as described in Sect. 12.31 ) or the B2C method, which 
is similar to the TLTD method. 
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Table B.l. Stellar and disk properties for young stellar objects. Column 10: the accretion rates are derived from the Her em ission line 
(see Sect. l3.2.Tb . besides objects SXCha, WXCha, and XX Cha, whose accretion rates are from Hartmann et al. ( 1998). Columns 
13, 14, 15: the average sizes of amorphous ((a)^ sil ) and crystalline grains ((a)^ sLsi] ) and the mass fractions of crystalline grains 

(f„ yst ) are derived from fitting silicate features around 10/mi with the TLTD or B2C methods. Columns 16, 17, 18: the average sizes 
of amorphous «a)^ m sil ) and crystalline grains ((a)^ sil ) and the mass fractions of crystalline grains (f^-y St ) <ire derived from fitting 
silicate features around 20-30//m with the TLTD or B2C methods. 
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Table A.l. Fractions of YSO with excess emission at wave- 
lengths less than ~8/im in different SFRs. Column 1: " Sparse 
stellar associa tions. Column 2: th e median age of each SFR es- 
timated from Baraffe et al. (1998). Column 3: the references for 
the median age of SFR. Column 4: the disk frequencies for all 
known members. Column 5: the references for disk frequencies 
of all known members. Column 6: the disk frequencies for the 
members with masses of O.1-O.6M . Column 7: the references 
for disk frequencies listed in Column 6. 
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